Mesoscopic fluctuations of the Coulomb drag at u=l/2 
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We consider mesoscopic fluctuations of Coulomb drag transresistivity between two layers at a 
Landau level filling factor v = 1/2 each. We find that at low temperature sample to sample 
fluctuations exceed both the ensemble average and the corresponding fluctuations at B = 0. At the 
experimentally relevant temperatures, the variance of the transresistivity is proportional to T -1 ^ 2 . 
We find the dependence of this variance on density and magnetic field to reflect the attachment of 
two flux quanta to each electron. 

PACS numbers: 73.43.Cd, 73.23.-b, 71.10.Pm 



Measurements of Coulomb drag between two parallel 
electronic systems at close proximity are a useful 

tool for studies of electron-electron interactions. In these 
measurements a current I\ is driven in one (" active" ) 
of the systems. As a consequence of inter-layer electron- 
electron scattering, a momentum q is transferred from the 
active system to the other (" passive" ) one with a voltage 
V-i developing in the latter. The ratio pr> — —V2/I1 is 
known as the transresistivity or the drag coefficient. 

For clean systems at zero magnetic field |6| -|8| p d °c T 2 
(with T being the temperature), as follows from Lan- 
dau's phase space argument. Recently, Coulomb drag 
was studied experimentally |^,|^] in a system of two two- 
dimensional electron layers in a strong magnetic field 
B. When the Landau level filling factor in each layer is 
v = 1/2, the effect is enhanced by several orders of mag- 
nitude relative to B = 0, and follows a sub-quadratic 
temperature dependence. Moreover, at very low tem- 
perature it was found to saturate to a non-zero value. 
Theoretically |@, Pd oc T 4 / 3 with a prefactor much 
larger than at B = 0. In the presence of disorder §|l|, 
both at B = and at v = 1/2, p D oc T 2 InT. 

Mesoscopic fluctuations of Coulomb drag at B = 
were recently studied theoretically [fr2| . When kpd 3> 1 
(hp is the Fermi momentum and d the inter-layer sep- 



aration), Coulomb drag is proportional to 



, dn v 



where 



a is the single layer conductivity and n is the electron 
density. On average, the derivative is a/n. Its fluctua- 
tions for a fully coherent sample are e 2 /TlEtN(Q) with 
E T = ^ being the Thouless energy and N(0) being the 
thermodynamic compressibility of electrons. Hereafter 
D is the diffusion constant of electrons. In the absence 
of phase breaking, for large enough samples the fluctua- 
tions exceed the average: (Pd)/{pd) 2 — (L/l) 4 , where I 
is the elastic mean free path. As usual, mesoscopic fluc- 
tuations of the drag result from quantum interference 
between different trajectories and thus are suppressed by 
phase breaking processes, leading to a non-trivial tem- 



perature dependence. 

In this paper we study fluctuations of Coulomb drag 
in the v = 1/2 case, where each layer is described E3] 
as a Fermi liquid of composite fermions. Each composite 
fermion is an electron with two attached flux quanta juj , 
that interacts with the others both electrostatically and 
by means of a Chern-Simons interaction. While most 
properties of the v = 1/2 state can be analyzed in terms 
of semi-classical dynamics of composite fermions, this 
study requires a quantum mechanical analysis. 

We start with a qualitative discussion, then sketch the 
calculation whose details are to be published elsewhere, 
and conclude by suggesting experimental probes to dis- 
tinguish the effect of mesoscopic fluctuations from other 
phenomena. All along, we use a and p to denote the 
measurable electronic conductivity and resistivity. Com- 
posite fermion conductivity, resistivity, and diffusion con- 
stant are denoted by o c $ , p c f , and D c ? . For simplicity 
we consider only the longitudinal drag conductivity and 
resistivity, denoted by crjj and po- We note in pass- 
ing that the Hall drag resistivity, which vanishes when 
disorder-averaged, would fluctuate with the same corre- 
lation functions as the longitudinal one, pe- 
lt is convenient to calculate mesoscopic fluctuations 
for the drag conductivity, ctd, rather than for pr>. In 
fact, at v = 1/2, to the lowest non- vanishing order in the 
screened interlayer interaction || , 

pd ~ f^r) °d, (1) 

In the same approximation, an can be expressed in terms 
of the non-linear susceptibility of the system T (to be de- 
fined below) and the propagator T> of the screened inter- 
layer Coulomb interaction: 

1 f duj ( d hio \ n „ a „ „ 
°d = — I — f — coth — V^V^Tl*, (2) 
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Here, S is the area of the sample, numerical subscripts 
indicate spatial coordinates, and are implied to be in- 



1 



tegrated over. The non-linear susceptibility of electrons 
Tfj(uj) is a response function relating a voltage V{ri)e lult 
to a dc current it induces by the quadratic response: 



J a = / dri / dr 2 r? 2 (w)V(n)V(r 2 ), 



(3) 



with J being the induced dc current. From gauge invari- 
ance / drxTf^uj) = J dr 2 T^ 2 (uj) = 0. 

The retarded (advanced) propagators of the screened 
interlayer interaction are calculated by means 

of a standard random phase approximation (see, e.g. 
Rcf. In momentum representation, 



V R(A) 



[n«W]" 



2ne 2 dU R ( A ) ' 



U R(A) 



N(0)Dq 2 
^fitu + Dq 2 



(4) 



In the diffusive regime the non-linear susceptibility T can 
be obtained from Ohm law, j — &E — eDVn, where E is 
the electric field. Combined with the continuity equation, 
Ohm law yields for the linear response of the density to 
an applied field, 



1 iq a a af3 E l3 (q,uj) 
e —iuj + Dq 2 



(5) 



where (. . .} indicate disorder-averaged quantities. 

A non-linear response to the electric field results 
from the density dependence of the conductivity J^ c = 
Re {^) n ( t l^ LJ ) E (-1^- LJ )^ and yields, 



/ r -n = 1 dja^) ^ uj<r xx q 2 



dn 



D 2 q 4 



e^Vanll* (6) 
on 



The density dependence of the conductivity is a measure 
of electron-hole assymetry, which is essential for Coulomb 
drag The fluctuations in pn result from mesoscopic 
fluctuations in §2. 

cm 

In the absence of a magnetic field T 1 is parallel to 
q. In contrast, at v = 1/2 the large Hall component 
of the conductivity leads to T 7 which is approximately 
perpendicular to q. In both cases the disorder-averaged 
conductivity is linear in the density, and d g n ks 
Substituting this approximation in Eq. (|^), and using 
Eqs. (0) and (0), we arrive at the familiar result B, 



(PD) 



2-irh 



3e 2 (k F d) 4 



hi 



(7) 



with T = 4na xx /ed. Note that Eq. (Q) holds both at 
zero magnetic field and at v = 1/2. Quantitatively, it 
yields very different results in the two cases, since the 
electronic longitudinal conductivity o~ xx at v = 1/2 is 
about three orders of magnitude smaller than at B = 0. 

We now turn to the discussion of mesoscopic fluctua- 
tions of pd- The first step is the estimate of fluctuations 



in r(q, us). At v = 1/2 the main source of these fluc- 
tuations is the derivative The rest of the parame- 
ters, i.e., the compressibility and diffusion constant, can 
be approximated by their average values (as in Eq. (||)), 
since their fluctuations are much smaller than their av- 
erage. To estimate the fluctuations of |^ we express this 
derivative in terms of response functions for composite 
fcrmions, which are the Fermi liquid quasi-particles of 
the v — 1/2 state. The conductivity matrices of elec- 
trons and composite fermions are related by 



2h 



(8) 



with e being the two dimensional anti-symmetric tensor. 

On average a c * , the composite fermion conductivity 
matrix at v = 1/2, is diagonal with both diagonal ele- 
ments being \g c f (where g c f is the dimensionless con- 
ductance of composite fermions) . In the limit of g c f ^ 1 , 
the electronic longitudinal conductivity is inversely pro- 

portional to g c f, a xx « xigj' and so ^ s a ^ so ^ ae e ^ ec ~ 
tronic diffusion constant. 

For a particular disorder realization the system is not 
isotropic, and a c f is not diagonal. In the limit of g c f ^> 1, 
Eq. (@) leads to 



8a 



Sa cf 



(9) 



where 5o~, 5o~ c f are the deviations of a, a c f of a particular 
sample from their average value. The fluctuations of ^ 
then stem from fluctuations of da c f /dn. This derivative 
is taken with the external magnetic field B kept con- 
stant. At the same time, the magnetic field experienced 
by the composite fermions, AB = B — 2<&on, is not con- 
stant: when the density is varied, AB varies as well. The 
derivative must then be taken with respect to both n and 
AS, 



da c f 



dn 



On 



2$ 



da c f 



AB 



dAB 



(10) 



We now estimate the fluctuations of T for a phase co- 
herent sample of size L. The scale for the variance of the 
conductivity a CJ at T = is — . The typical magnetic 
field scale is AB ~ $q/L 2 , therefore the estimate for the 



fluctuations of 2<E>o 



da cf 
dAB 



is 4-L 2 . 

h 



The fluctuations of 



da af 
dn 



are smaller by a factor of g c f, and can there- 
fore be neglected. Consequently, the variance of both 
components of T is 



ST ~ iq- 



ImlT 



(11) 



L 

9c f, 

2 

b = the fluctuations of T 

(^r) 4 ' which in the diffu " 
sive regime is much larger than unity. This by itself is a 



d{c 



dn 



Since on average 



can be re-written as |p^ 2 



2 



measurable conclusion: In the diffusive regime in a fully 
coherent sample the fluctuations of the acousto- electric 
current are much larger than its average. Equation ( |TT| ) 
holds as long as the thermal length Lyf = ^h,D c f /T and 
the phase breaking length L C J are much larger than L. 

Having now an estimate for the fluctuations of T, we 
estimate the fluctuations of po as: 



Spi 



T 



dtv 



dq ST 2 \V(q,u)\ 2 (12) 



| 9 |>l/min(L,L c /) 



As commonly happens in mesoscopic fluctuations, the in- 
tegral over q is dominated by its lower limit. We start 
with the case where the phase breaking length is the 
largest scale in the problem. We describe the temper- 
ature dependence of the transresistance, and denote the 
r.m.s of the fluctuations by Sp^. 

At extremely low temperatures T <C (nd)hD / L? , all 
the excitations with energies smaller than T contribute 
and we obtain with the help of Eqs. (|l2|) and (0) 



where the coefficients a 1.2 are of order unity and we were 
able to calculate a 3 w 0.2(32/9tt) = 0.23. 

While the actual magnitude of the mesoscopic fluc- 
tuations depends on the precise source of phase break- 
ing, their temperature dependence is robust. All generic 
models of phase breaking in two dimensions lead to 
1/t v oc T, so that the temperature dependence follow- 
ing from Eq. (0) is (5p 2 D ) ocl/T. 

The dependence of the mesoscopic fluctuations on g c f 
requires further specification of the model. Our prelim- 
inary study suggets that the main mechanism for the 
phase breaking is the quasi-elastic scattering of compos- 
ite fermions from the thermal quasi-static fluctuations of 
the Chern-Simons magnetic field, (2<& )~ 2 (-B(r)i?(r')) = 



(Sn(r)Sn(r 1 )) « S(r 



-n tn(0) _ 

K,d 



e being the bulk dielectric constant). Due to the accu- 
mulation of the Aharonov - Bohm phase from such fluc- 
tuations, (Lpf) 2 is the area at which the thermal fluc- 
tuations of the electron number are of the order of one 
ii \<L cf drdr' (n(r)n(r')) ~ 1, which result in 



S P °n 



1 



e 2 M 2 




(nd)El 



•cf 
T 



9lf 



(13a) 



where El 



cf 



h~D c f/L 2 is the Thouless energy for 
the composite fermions. At higher temperature 
(Kd)hD/L 2 <T< E%f, the processes with the energy 
transfer ui > (nd)D/L 2 are suppressed, leading to 



5 P ° L 



T 



D 



e 2 {nd)g 2 f El 



cf 



{hd)E T f <T<E cf ( 13b ) 



At yet higher temperature, T > Ej! , ST itself is tem- 



perature dependent, and it is suppressed by a factor 
of y E!p /T similar to conductance fluctuations, see e.g. 



L6| . This yields a temperature independent result 



SP° D 



1 



M)3c/ ' 



T > E C J . 



(13c) 



Let us now discuss the effect of a finite phase breaking 
length L C J . As is well known for mesoscopic fluctuations, 
averaging in a large sample, L 3> L C J , is carried out by 
summing over the statistically independent contributions 
of (L/L C J) 2 patches of a size L C J each. The contribution 
of each patch is given by Eqs. ([l3]), where the Thouless 
energy E^ is set to h/r^f = hD c f/(L c J) 2 . As a result of 
the averaging {5p 2 D ) = (S p° D ) 2 {L C J / L) 2 . Using Eqs. 
we obtain 



(P 



=4 n 4 




1, OL\ 



ndh 



x mm 



a 3 ,a 2 (TT° f /hy 



(14) 



2iTe 2 d 



cf 



9cfT 
nd 



(15) 



This estimate can be understood very simply: a thermal 
fluctuation of a charge e in one layer and — e in the other 
layer over a scale L C J involves a charging energy cost of 

TTWvi ■ The energy available for that fluctuation is T. 
Balancing the two energies determines L C J . 




Substituting estimate (15) into Eq. (|14| ) and using the 
condition g c ? 3> nd, we find up to a numerical constant 



(16) 



Our estimate ( p"5| ) for the phase breaking rate implies 
Tt v <^ h. This by no means indicates a collapse of the 
Fermi liquid picture of composite fermions, since most of 
the phase breaking results from scattering off the Chern- 
Simons field fluctuations whose dynamics [with charac- 
teristic frequency T/(h~g c f)] is very slow compared to t v , 
but fast compared to the time of the experiment. Field 
fluctuations which are static on the scale of the experi- 
ment time affect the mesoscopic fluctuations only by af- 
fecting g c f. Field fluctuations that are faster than that 
scale make the potential landscape seen by the composite 
fermions time dependent, and lead to a suppression of the 
mesoscopic fluctuations by partial ensemble averaging. 

We now use Eq. ([l6]) to estimate the value of the fluc- 
tuations of the transresitance of a realistic sample. At 
T = 0.6K the average drag Q (p D ) = 15fi/D with the 
interlayer spacing d = 300A For that sample, the single 
layer resistance at v = 1/2 is B R = 3fcS!/n, which gives 



3 



the value g c f ~ 8. Given the sample's size L ~ 100/ito, 
we estimate the magnitude of the fluctuations Eq. Jl6| ) 
as Spn ~ 0.3f2. For lower temperatures and smaller sam- 
ples we expect the fluctuations to be stronger, possibly 
exceeding the average. Increasing current supresses the 
mesoscopic fluctuations mostly by heating of the elec- 
trons. 



J a 




FIG. 1. Electronic non-linear response within the compos- 
ite fermion RPA. Wavy lines are gauge field propagators Jl3{ . 
Bubbles are composite fermions response functions. 

We now turn to sketch the calculation, whose de- 
tails will be published separately. Within the composite 
fermion random phase approximation (RPA), T is de- 
picted diagrammatically on Fig. [TJ. The center triangle 
is the non-linear response of the composite fermions T c ' 
to a driving field. Since composite fermions interact with 
scalar and vector potentials, r c ^ is a tensor with compo- 
nent (r c ^)^ /3 , where a/3 are the directions of the driving 
fields, and p is the direction of the induced current. 

In the limit of g c f 3> 1 the analytic expression corre- 
sponding to the diagram on Fig. [TJ is approximated by 



e 2 g cf 1 J34 ^ V 



where numerical subscripts indicate spatial coordinates, 
we imply integration over ^4, and the polarization oper- 
ators are defined in Eq. ([I]) . The mesoscopic fluctuations 
in r result from fluctuations in r c /. Within the RPA, 
the latter is approximated by the corresponing response 
function for non-interacting particles: 

FcfT^ = J^ [J&., M uM + J2W(-w,e)] . (17) 

J?o-nn(wi e ) = I tanh - — — — tanh — 

x [Gf 2 (e - u>) - G? 2 (e - u )] 3, [G R (e)j«G A (e)] 21 j u ; 

The further analysis of the statistics of the vertex (|l 
follows the lines of Ref. p2| , and results in Eq. (fj~ 



Closing the paper, we point out that two analy- 
ses of mesoscopic fluctuations of the Coulomb drag 
near v = 1/2 are of interest. First, one varies the 
magnetic field and measures the correlation function 
(p D {B)p D (B + SB)) - (p D (B))(p D (B + SB)), with B 
close to the v = 1/2 value. An experimental study of 
the decay of this correlation function is a way to mea- 
sure L C J: the characteristic magnetic field of the decay is 
SB* ~ $ /[££/] 2 ■ Second, one varies the electron density 
in one of the layers and measures the correlation function 
(p D (n)p D (n + Sn))-(p D (n))(p D (n + Sn)). Again, the de- 
cay of this function is governed by L^f: the characteristic 
density change 5n* at which it decays is expected to cor- 
respond to half of an electron in a phase coherent region, 
i.e. Sn* = 1/2[L£/] 2 . This statement should hold as long 
as the composite fermion cyclotron radius is much larger 
than its mean free path, i.e., for \v — 1/2 < (2g c f)^ 1 . 

It is noteworthy that near zero magnetic field SB* = 
&o/L^, but Sn* — kfl/Lp. Since at B = the electrons 
do not carry flux tubes, the density change Sn* has to be 
such that the chemical potential changes by Ti/t v . Thus 
an experimental observation of = 2<5> would in some 
sense be a verification of the attachment of flux to the 
Fermi liquid quasi-particles at and close to v = 1/2. 
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